function [act logact] = gmmactivRobustUncertainty(mix, data)


K = mix.ncentres;
N = size(data.Y, 1);

act = zeros(N, K);
logact = zeros(N, K);
for k = 1:mix.ncentres
    [act(:, k) logact(:, k)] = gaussRobustUncertainty(mix.centres(k, :), mix.covars(:, :, k), data.Y, data.C);
end
